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How may one locate a black hole? This question may 
be quite challenging even for numerical computation, for 
the location of a black hole horizon is a delicate issue. 
Here we propose a method which generically gives the 
precise location of any stationary horizon and should give 
the approximate location of a nearly stationary horizon. 
Hence it should be useful for numerical computations in 
general relativity. 

For a general black hole spacetime, the location of its 
surface (the event horizon, which is the boundary of the 
region from which causal curves can go to asymptotic fu¬ 
ture null infinity) depends on the future evolution of the 
spacetime and is not determined locally. However, for 
black holes that settle down to become stationary, one 
might ask whether one can find a local invariant that is 
generically nonzero off the horizon but vanishes on the 
horizon. In particular, one might look for a scalar polyno¬ 
mial curvature invariant , which is a scalar obtained 
by complete contraction of all the indices of a polyno¬ 
mial in the Riemann curvature tensor Rap-yS and its co¬ 
variant derivatives. For example, Anders Karlhede, Ulf 
Lindstrom and Jan Aman Q showed that R°‘^^^’'^Rap-ys-,e 
crosses zero and switches sign as one crosses the horizon 
of the Schwarzschild metric, and one can easily show that 
this is also true for any smooth static spherically symmet¬ 
ric horizon. An invariant that vanishes on more general 
stationary horizons could be useful in numerical relativity 
for finding the approximate location of the horizon once 
the spacetime has settled down to become approximately 
stationary. 

Majd Abdelqader and Kayll Lake Q have recently 
found a local scalar polynomial curvature invariant that 
vanishes on the horizon of the Kerr black hole. After 
casting this invariant into a simpler form that is propor¬ 
tional to the squared norm of the wedge product of two 
curvature-invariant gradients, we realized that the proce¬ 
dure generalizes to give a way to locate any nonsingular 
stationary horizon in terms of local curvature invariants. 
Essentially, if one constructs as many gradients of inde¬ 
pendent curvature invariants as the local cohomogeneity 
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of a generic stationary spacetime, at a generic point in 
the spacetime outside an horizon these gradients will be 
linearly independent and spacelike, but on the horizon, 
a linear combination will become null. This implies that 
the squared norm of the wedge product of the gradients 
will vanish on a stationary horizon. 

Our procedure does not assume the Einstein equations, 
but only a smooth stationary horizon (a null hypersurface 
generated by a Killing vector held that is timelike out¬ 
side the horizon but null on the horizon). Therefore, in 
principle the invariants to be used in the procedure could 
include those from the Ricci tensor R^p , such as the Ricci 
scalar R = R’^ or R°'^Rap. However, because the main 
interest may be in Ricci-hat (Rap = 0) spacetimes that 
solve the vacuum Einstein equations, in the examples be¬ 
low we shall consider only invariants obtained from the 
Riemann curvature tensor Rap-yS (or from the Weyl ten¬ 
sor Coip-yS, which is the trace-free part of the Riemann 
tensor) that need not vanish even when R^p = 0. 

In the original (30 Dec. 2014) version of their paper, 
Abdelqader and Lake @ gave the following six curva¬ 
ture invariants for the Kerr metric, which we here copy 
directly from that paper: 
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where Cap~fS is the Weyl tensor, C*apyS its dual, = 
Ii , and A . Then they showed that in 

the Kerr metric the dimensionless scalar nonpolynomial 
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curvature invariant (with scalar polynomial curvature in¬ 
variant numerator) 

^ _ (/5 + hf - (12/5)2 (/i2 + (732 + 7^2) 

^2 = - ^ - UJ 

108 [h’^ + h'^f 

vanishes on the black hole horizon. 

We found that the syzygy (in this context a functional 
relationship between invariants for a particular spacetime 
geometry) that Abdelqader and Lake @ discovered for 
the Kerr metric, that 

/6-/5 + y(/l/3-/2/4)=0, (8) 

may be expressed as the real part of the complex syzygy 
(1^1 + *^ 2 ) (/i -I- i/ 2 ) = — (/i + i/ 2 ) (-I 3 + Hi) 

( 9 ) 

whose imaginary part gives a syzygy for a seventh invari¬ 
ant, 

/r = fc^Z^ = ^(/i/4 + /2/3) • ( 10 ) 

Using both the real and imaginary parts of this complex 
syzygy, one may readily see that 

27 (/i2+ Q 2 = {k^kf^){ijn - {k^mi.ki , 

(11) 

which is the squared norm of the wedge product dli A dl 2 
of the gradients of the Kretschmann invariant Ii and of 
the Chern-Pontryagin invariant I 2 (for a Ricci-flat space- 
time; otherwise in four dimensions the Kretschmann in¬ 
variant is = /i+ 2 /?a/ 3 /?“^ —(l/ 3 )i? 2 ^ though 

the Chern-Pontryagin invariant R*a/ 3 -yS = I 2 has 

no correction from the Ricci tensor). 

Although it is incidental to the main point of our pa¬ 
per, we also found another syzygy from a minor extension 
of the work of Abdelqader and Lake Q, 


(/l + z/2)" (/3 - ihf 
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where S denotes the imaginary part of what follows. 
(The real part is equal to the absolute value in this case 
and does not vanish to give yet another syzygy.) With 
these three syzygies for the seven real scalar polynomial 
curvature invariants, /i,..., /y, one is left with four inde¬ 
pendent invariants, precisely the right number to deter¬ 
mine generically (up to eight sign choices) the values of 
the two Kerr parameters M and a and of the two nonig- 
norable coordinates r and 9. (None of these three remain 
syzygies when one adds a cosmological constant, so the 
syzygies for that case remain to be found.) 

Eq. (fTTl) for the Abdelqader-Lake invariant Q 2 that 
they found vanishes on the horizon of the Kerr metric in¬ 
spired the realization that the squared norm of the wedge 
product of n gradients of independent local smooth cur¬ 
vature invariants would vanish on the horizon of any sta¬ 
tionary black hole, where n is the local cohomogeneity 


of the spacetime (the codimension of the maximal di¬ 
mensional orbits of the isometry group of the local met¬ 
ric, ignoring the breaking of any of these local isome¬ 
tries by global considerations, such as the way that the 
identifications of flat space to make a torus break the 
local rotational part of the isometry group), so that at 
generic points ontside an horizon of a generic metric of 
local cohomogeneity n, the squared norm of the wedge 
product would be positive. (We use the metric signature 
—!-+•••.) Let us now prove that this squared norm 
vanishes on a stationary horizon. 

Let D be the dimensionality of the spacetime, and let 
m be the maximal dimension of the orbits of the local 
isometry group. Then n = D — m is the local cohomo¬ 
geneity of the spacetime. For a stationary spacetime of 
local cohomogeneity n with an event horizon, let de¬ 
note the Killing vector field that on the horizon is its 
null generator, which is orthogonal to the null horizon 
hypersurface. In a neighborhood outside the horizon, 
will be timelike, with < 0 in our choice of signature. 
Let for i ranging from 1 to n, denote a set of n 

functionally independent nonconstant scalar polynomial 
curvature invariants (total scalar contractions over all in¬ 
dices of polynomials of curvature tensors and of their co¬ 
variant derivatives). For example, some or all of the 
could be chosen from the set R.. .Ij given above (though 
remembering that I 2 , R, h, and R above are only de¬ 
fined for /) = 4), but one could also choose invariants 
from total scalar contractions of other polynomials in the 
Riemann curvature tensor and its covariant derivatives. 

Then let f,e the set of exterior derivatives (gra¬ 

dients) of the n scalar polynomial curvature invariants, 
with components Since the curvature 

invariants are invariant under translations by the lo¬ 
cal isometry group, their gradients all lie with the n- 
dimensional local cohomogeneity part of the cotangent 
space at each point The wedge product of n gradi¬ 
ents, the n-form W = dS^^'> A dS^^'> A • • • A dS^'^'>, will be 
proportional to the volume form in this n-dimensional 
part of the cotangent space at each point of spacetime, 
with a proportionality depending upon the spacetime 
point. At generic points in a generic spacetime, the pro¬ 
portionality will be nonzero where the gradients of the n 
scalar invariants are linearly independent, but there can 
be a set of points (generically hypersurfaces) where the n 
gradients are not linearly independent, so that the wedge 
product vanishes there and the proportionality is zero. 

The Hodge dual *W is then an m-form that at 
each point of the spacetime lies entirely within the m- 
dimensional part of the cotangent space that is gener¬ 
ated by the local isometries. It will be proportional to 
the m-dimensional volume element of that part of the 
cotangent space, though the location-dependent propor¬ 
tionality will be zero at the same set of spacetime points 
at which the wedge product W vanishes. 

On a stationary horizon, which is a null hypersurface 
generated by the Killing vector that is timelike outside 
the horizon but null on the horizon (and hence both par- 
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allel and orthogonal to the horizon), the m-dimensional 
part of the cotangent space will include the null one- 
form, with one-form components , which is 

metrically equivalent to the Killing vector that is null 
and hypersurface-orthogonal at the horizon, as well as 
other spacelike one-forms metrically equivalent to other 
Killing vectors at the horizon that are all spacelike if 
TO > 1. Therefore, this m-dimensional part of the cotan¬ 
gent space will be null on the stationary horizon, and 
hence *W will also be null, having zero squared norm. 
Since the squared norm of W, appropriately defined, is 
the same as that of its Hodge dual *W, the squared norm 
of the wedge product W oin scalar polynomial curvature 
invariants will be zero on a stationary horizon. 

There is of course in this argument the implicit as¬ 
sumption that the spacetime is sufficiently regular at the 
event horizon that the gradients of the n scalar polyno¬ 
mial curvature invariants are well-defined at the horizon, 
so that their wedge product W is well-defined and finite. 

We summarize this result by the following theorem: 

Theorem. For a spacetime of local cohomogeneity n that 
contains a stationary horizon (a null hypersurface that is 
orthogonal to a Killing vector field that is null there and 
hence lies within the hypersurface and is its null gener¬ 
ator) and which has n scalar polynomial curvature in¬ 
variants whose gradients are well-defined there, the 
n-form wedge product W = dS^^^ A dS^'^'^ A • • • A has 

zero sguared norm on the horizon, 

( 13 ) 

where the permutation tensor is -|-1 if ai... an 

is an even permutation of Pi... Pn, is —1 if ai ... an is 
an odd permutation of Pi... Pn, and is zero otherwise 
(including all cases of repeated indices upstairs or down¬ 
stairs). 

Then ||W|p is itself a scalar polynomial curvature in¬ 
variant that vanishes on any stationary horizon smooth 
enough for W to be well-defined there. Of course, for this 
to be useful for locating a stationary horizon, the n scalar 
polynomial curvature invariants should be chosen to 
be functionally independent so that ||lT|p is positive at 
generic points in the spacetime outside the horizon. This 
will not always be possible, such as for the cosmological 
horizon of de Sitter spacetime, which is totally homoge¬ 
neous and hence has cohomogeneity n = 0 , though of 
course for this spacetime there are cosmological horizons 
running through every point. 

This raises the following question: For a spacetime of 
local cohomogeneity n with a Killing horizon that is a 
locally unique hypersurfaccQ: are there always enough 


r By locally unique, we mean that any neighborhood of the horizon 


functionally independent scalar polynomial curvature in¬ 
variants so that one can choose n of them in such a way 
that the squared norm of the wedge product of their gra¬ 
dients is generically nonvanishing away from the horizon? 
If the answer is ‘no,’ then although our theorem above 
would remain true, the squared norm of the wedge prod¬ 
uct of the gradients of any n scalar polynomial curvature 
invariants would vanish in a whole neighborhood of a sta¬ 
tionary horizon, so it would not be useful for locating the 
horizon. Based upon the work of Alan Coley, Sigbjprn 
Hervik, and Nicos Pelavas one might conjecture 

that the answer to the question is ‘yes,’ but so far we do 
not have a rigorous proof. 

Let us now consider various examples in four¬ 
dimensional spacetimes of the fact that the squared norm 
of the wedge product of n gradients of curvature invari¬ 
ants vanishes at the horizon of a stationary black hole. 

First, consider a spherically symmetric static black 
hole, so the codimension is n = 1 (the radial direction, 
say with coordinate r). Then any smooth curvature in¬ 
variant, such as the Kretschmann invariant Rap-ysR°‘^^^, 
which is the same as Ii for a Ricci-flat spacetime, will 
have a gradient in the radial direction that becomes null 
on the horizon (e.g., as seen in a frame parallely propa¬ 
gated by an infalling geodesic), so the square of the norm 
of the gradient of the Kretschmann invariant (with only 
one gradient in the wedge product in this n = 1 case) 
will be zero on any spherically symmetric static black 
hole horizon. 

Second, consider a stationary axisymmetric black hole 
in four dimensions, which has two commuting Killing vec¬ 
tor fields and cohomogeneity n = 2, such as the Kerr 
metric. Then we need the wedge product of two gradi¬ 
ents for its squared norm to be positive at generic loca¬ 
tions outside the black hole but vanishing on the hori¬ 
zon. For a rotating black hole such as Kerr, one can 
take the squared norm of the wedge product dli A d /2 
as Abdelqader and Lake @ effectively did (though origi¬ 
nally without realizing explicitly that their invariant Q 2 
is proportional to the squared norm of this wedge prod¬ 
uct before we discovered and pointed out this fact). For 
a static axisymmetric black hole, the invariants I 2 , I 4 , 
Iq, and Ir vanish, so one could instead take the squared 
norm of d/i Ad/ 3 , of d/i Ad/ 5 , of dl^Adl^, or of any other 
pair of independent nonvanishing scalar polynomial cur¬ 
vature invariants, such as the trace of higher powers of 
the curvature tensor like 

In these cases with more than one Killing vector that 
all commute (so that the dimensionality to of the isome¬ 
try group is the same as the number of Killing vectors), 
the squared norm of the wedge product of n independent 
scalar polynomial curvature invariants will also vanish 


includes points not on a Killing horizon, unlike the case for de 
Sitter or certain other examples, such as spacetimes having a 
covariantly constant null vector field, that have Killing horizons 
through all points. 
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at the fixed points of the Killing vectors other than the 
one that becomes null on the horizon, such as on the 
axes of axisymmetric spacetimes. In such cases we do 
not know how to construct scalar polynomial curvature 
invariants that vanish only on the horizon, so that will 
be left as a challenge for the future, as well as the chal¬ 
lenge as to whether for specific spacetimes such as Kerr, 
there is a scalar polynomial curvature invariant without 
any covariant derivatives of the Riemann curvature ten¬ 
sor that vanishes on the horizon or that vanishes only on 
the horizon. 

Third, consider a distorted static black hole in four 
dimensions that has no spatial Killing vector fields on 
and outside the horizon, so the local cohomogeneity is 
n = 3. Then we need the squared norm of the wedge 
product of three gradients of scalar polynomial curvature 
invariants, such as dli A dl^ A dl^ or wedge products 
using the gradients of the traces of higher powers of the 
curvature tensor and/or of its covariant derivatives. 

There are of course analogous examples in spacetimes 
of higher dimension D > A. If they are static and spheri¬ 
cally symmetric, then for any dimension H > 4, the local 
cohomogeneity is n = 1 , and the gradient of any scalar 
polynomial curvature invariant (such as the Kretschmann 
invariant) that has a nonvanishing gradient just outside 
the horizon will have its squared norm becoming zero on 
the horizon. 

The general Kerr-NUT-(A)dS metric in spacetime di¬ 
mension D 0 has local cohomogeneity n = [D/2\, the 
greatest integer not greater than D/2. That is, the D = A 
case has two commuting Killing vector fields and local co¬ 
homogeneity n = 2 (as the special case of the Kerr metric 
does with both the NUT parameter and the cosmological 
constant zero). Analogously, the D = 5 case has three 
commuting Killing vector fields (one becoming null on 
the horizon and two axial ones) and also n = [5/2J = 2 . 
Then D = 6 and D = 7 have n = 3, D = 8 and D = 9 
have n = A, D = 10 and D = 11 have n = 5, etc. In 
each case the squared norm of the wedge product of the 
gradients of n scalar polynomial curvature invariants will 
vanish on the horizons (as well as on the axes surfaces, 
the fixed points of the axial Killing vector fields). 

One might conjecture that one could take the traces 
of the 2nd through the (n -I- l)th powers of the Riemann 
curvature tensor as one choice of n invariants. We have 
looked at sample values of the 3 parameters (mass M, 
rotation parameter a, and NUT parameter L at fixed 
cosmological constant) of the Kerr-NUT-(A)dS metric 
in spacetime dimension D = A and have found that in¬ 
deed the square of the norm of the wedge product of the 
traces of the 2nd and 3rd powers of the Riemann curva¬ 
ture tensor (the Kretschmann invariant and 

is generically nonzero away from the 
horizons and axes, but there are curves in the (r, 0) plane 
(hypersurfaces of the full spacetime that are neither hori¬ 
zons nor axes) where the square of the norm of the wedge 
product vanishes. 

In particular, when the NUT parameter L is zero, the 


Kerr-(A)dS spacetime has a reflection symmetry about 
the equatorial plane hypersurface, if one ignores the 
change in sign of the volume element given by the to¬ 
tally antisymmetric Levi-Civita tensor. In this case, gra¬ 
dients of scalar polynomial curvature invariants that do 
not involve odd powers of the Levi-Civita tensor (e.g., not 
counting / 2 , / 4 , and m. D = A) will lie in the (n — I)- 
dimensional local cohomogeneity part of the hypersurface 
of reflection symmetry, so that the wedge product of any 
n such gradients will vanish there. When the NUT pa¬ 
rameter is nonzero, there is no reflection symmetry, but 
generically there remain distorted hypersurfaces where 
the wedge product vanishes. 

We found that there are generically other hypersur¬ 
faces away from the horizons and axes on which the 
wedge product of n = 2 gradients vanishes. How¬ 
ever, we found that if we took the three scalar polyno¬ 
mial curvature invariants = 

and 5(3) = then we 

have three pairs of wedge products, IU 12 = dS^^^ A dS^‘^\ 
IU 23 = d5(2) A d5(3), and IU 31 = dS<^^'> A dS<-^\ that 
each vanish on certain curves in the (r, 9) plane and with 
each pair concurrently vanishing at certain intersection 
points in the (r, 6) plane, but with no triple intersec¬ 
tions where all three wedge products concurrently vanish 
for the sample values of M, a, and L that we chose. 
Therefore, it appears that if one takes the sum of the 
three squared norms of the wedge products each multi¬ 
plied by positive coefficients, then for generic parameters 
this will be nonzero everywhere except at the horizons 
and fixed points of the isometry. However, we did find 
that for certain choices of the Kerr-NUT-(A)dS parame¬ 
ters, even this sum of squares vanishes at certain sets of 
points in the spacetime away from the stationary hori¬ 
zons and fixed points of the isometries. We believe this 
vanishing anywhere could be eliminated for any set of the 
three parameters by choosing a sum of six squared norms 
of the wedge products of the six pairs of a suitable set of 
four gradients of scalar polynomial curvature invariants. 

One general procedure for calculating a scalar polyno¬ 
mial curvature invariant /(p), vanishing on a stationary 
horizon, that should generically be positive wherever the 
Killing vectors span the maximal dimension m of the lo¬ 
cal isometry group that includes a timelike direction (so 
that one is not at a stationary horizon or fixed point of 
the isometry, and so that the n independent directions 
of the local cohomogeneity are all spacelike) would be 
to take a set ol n + p scalar polynomial curvature in¬ 
variants, {5(*)} for i running from 1 to n -|- p with suf¬ 
ficiently large p. (For example, p could be the number 
of parameters for a certain class of spacetimes, such as 
p = D — 1 for the Kerr-NUT-(A)dS family of metrics [l^ 
of even dimension U, or p = £) — 2 for odd D, for which 
there are n + p algebraically independent scalar polyno¬ 
mial curvature invariants.) Then there are (n-|-p)!/(n!p!) 
ways of choosing n gradients to form the squared norm 
of wedge product n-form, and one can multiply each by 
the product of the squared norms of the remaining gradi- 
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ents and then sum over the (n+p)!/(n!p!) choices to get 
an invariant /(p) that will generically be positive where 
all the gradients are spacelike, as they generically would 
be under the conditions above. Of course, precisely what 
scalar polynomial curvature invariant /(p) one gets would 
depend upon which set of n + p scalar polynomial cur¬ 
vature invariants one chooses, and some special 

choices could still conceivably give /(p) = 0 even away 
from stationary horizons and fixed points of the isome¬ 
try group. 

In conclusion, we have found a way to generalize the 
results of Abdelqader and Lake @ for the Kerr metric 
to give ways of constructing scalar polynomial curvature 
invariants that vanish on any stationary horizon. These 
invariants might be useful for numerically estimating the 
location of horizons once a spacetime settles down to be¬ 
coming nearly stationary. They may also be interest¬ 
ing in nonstationary black hole spacetimes for defining 
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